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Abstract 

We show that a type of linear superposition principle works for several nonlinear differential equa- 
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1 Introduction 



The fact that Jacobi ehiptic functions arise naturahy as travehng wave solutions of many nonhnear 
systems has been known for quite some time (see for example Q). Although for the solitary wave solu- 
tions of these nonlinear equations, there is no superposition principle (except when the solitary waves 
are far apart), for the periodic solutions the situation turns out to be quite different. It has recently 
been shown Q that certain specific linear combinations of known periodic solutions of the Korteweg-de 
Vries (KdV) and modified Korteweg-de Vries (mKdV) equations as well as X(p^ theory, also satify these 
equations. This unexpected result is a consequence of some remarkable, recently established, identities 
involving Jacobi elliptic functions [Q. Basically, the identities take the cross terms generated by the 
nonlinear terms in the differential equations and convert them into a manageable form. The purpose 
of this article is to show that such a procedure also works for other well-known nonlinear equations, 
namely the Kadomtsev-Petviashvili (KP) equation, the nonlinear Schrodinger equation(NLSE), the 
sine-Gordon equation and the Boussinesq equation. It should be noted that the above list includes 
both integrable as well as nonintegrable systems. These equations are of interest in several diverse 
areas of physics. The NLSE governs the propagation of an electromagnetic wave in a glass fiber, or 
the spatial evolution of an electromagnetic field in a planar waveguide. Temporal solitons described 
by the NLSE were first observed in 1980 [^, and the first confirmation and studies of spatial solitons 
in planar waveguides were reported in 1988 ^. Similarly, the and the sine-Gordon equations 
arise in several condensed matter physics applications. 



2 The Kadomtsev-Petviashvili (KP) Equation 

The KP equation is a two-dimensional generalization of the KdV equation and is given by 

(Ut - 6UUx + Uxxx)x + ^Uyy = . (l) 

Properties of the KP equation are discussed in many texts Q. In particular, the simplest, periodic, 
cnoidal traveling wave solution is 

ni(x, y, t) = — 2a^dn^(^i, m) -|- /?a^ , = a{x + 'jay — bia^t) , (2) 

where a, 7, m and (3 are constants, and the "velocity" 61 is given by 

61 = 8 - 4m - 6/3 37^ . (3) 
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In this article, for Jacobi elliptic functions, we use the standard notation dn m), sn m), cn m), 
where m is the elliptic modulus parameter (0 < m < 1). The solution (2) remains unchanged when x is 
increased by 2K{m)/a, where K{m) is the complete elliptic integral of the first kind In the limiting 
case m = 1 (and /3 = 0), one recovers the familiar single soliton form — 2a^sech^(a(2; + 'jay — bia'^t)). 

We will make suitable linear combinations of solution (2) and show that the result is also a periodic 
solution of the KP equation. Our procedure consists of adding terms of the kind given in (2) but 
centered at p equally spaced points along the period 2K{m)/a, where p is any integer. The p-point 
solution is 

^ 2(i — l)K(m) 
Up{x, y, t) = —20? ^ + (ic? ; di = dn[,^p H , m] , = a{x + ^ay — bpO^t) . (4) 

i=i P 
Clearly, p = 1 is the original solution, but for any other p, we have new expressions which, as we shall 
show, also solve the KP equation. For convenience, we define the quantities Si and q in analogy to 
the quantity di defined above: 

2{i-l)K{m) , 2{i-l)K{m) ^ 

Si = sn^p + ^ ,m , a = cn[Cp + ^ \ m] . (5) 

p P 

The KP equation contains the KdV operator ut — 6uux + Uxxx- It has been shown in detail in ref. 
1^] that eq. ^ with 7 = is a solution of the KdV equation. The proof is based on the identity 

X: d^d] = Aiip, m) j2 dl + A2{p, m) . (6) 

i<j i=l 

This is one of many powerful new identities Q which reduce by 2 (or a larger even number) the degree 
of cyclic homogeneous polynomials in Jacobi elliptic functions. The constants Ai{p,m) and A2{p,m) 
in identity can be evaluated in general by choosing any specific convenient value of the argument 
^ of the Jacobi elliptic functions. The results for Ai{p,m) for small values of p are: 

-2(m -l + q'^^ 



Ai{p = 2,m)=0 , Ai{p = 3,m) = ^— , Ai{p = 4,m) = -2^1 - m , (7) 

1 — (7^ 

where 

q = dn(2K(m)/3,m) . (8) 
The limiting values at m = 0, 1 are also particularly simple: 

Ai{p,m = 0) = -^{p-l){p-2) ; Ai{p,m = 1) = . (9) 
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Figure 1: ui{x,y,t = 0) vs. for m = 0.95, a = 7 = 1, /3 = 0. 

Taking expression @ and using identity (^), the left side of the KP equation (|l|) becomes 

Ama^ {8 - 4m - 6f3 - b„ + 12 Ai(p,m)}—y^ SiCidi + 12mja^ — VsiQfii . (10) 

ax ay 
1=1 ^ 1=1 

Clearly, this vanishes if the velocity is given by 

6p = 8-4m- 6/3+ 12^1 (p, m) + 37^ . (11) 

Thus for this choice of velocity, the KP equation is solved by our p-point expression (Q). Effectively, our 
solutions of the KP equation and the corresponding solutions of the KdV equation have a difference of 
37^ in their velocities bp. Note that as in the KdV case, the results for bp can be positive or negative 
depending on the values of the parameters [^. The behavior of ui for the parameters m = 0.95, 
a = 7 = l, /? = Ois shown in fig. ||. 

In addition to the solution (^), another well known periodic solution of the KP eq. (|l|) of period 
AK{m) is 



vi{x,y,t) = cP- 



msn^ (r\\ , ra) it \pmcii{r]\ , m)dn(?7i , m) 



, ri\ = a{x ^ ^ay - q\o?t) (12) 



with velocity q\ = (—1 — m + 37^). 

Starting from this solution, we can again obtain many additional periodic solutions of the KP 
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equation of period 4K(m)/p in case p is an odd integer. The general p-point solution is given by 

p 

Vp{x,y,t) = a^^'^ [msj :t y/mcidi] , p odd, (13) 



i=l 

where we define 



4(i-l)K(m) ^ , ^ 4(i-l)K(m) ^ ~ ^ ^ 4(i - l)K(m) , 

Si = sn ??p H , m , a = cn H ,m\ , di = dn % H , m . 

p p p 

(14) 



As has been shown in detail in ref. |]2|, eq. (13) with 7 = is a solution of the KdV equation with 
velocity 

Qp = -(1 + m) - 6[Bi{p,m) - Ci{p,m)] , (15) 

where the quantities Bi{p,m) and Ci(p, m) come from the following identities: 

p p p 

m^SiSj = Bi{p,m) , m ^ SiSjSk = Ci{p,m) . (16) 

i<j i<j<k 1=1 

It is easily checked that even in the KP case (7 7^ 0), eq. (|l^) is an exact solution, the only difference 
being that the velocity in the KP case is larger by 87^. As an illustration, for the p = 3 case, it is 
easily shown that 

Bi(3,m) = -(l-g2) , C7i(3,m) = -m/(l-g2) , (17) 

so that the velocity of the KP soliton is given by 

q^ = -l-m + 6il-q^)--^ + 3j\ (18) 

1 — 

where q has been defined in eq. (^). Note that for p even, we do not obtain any new solutions. To 
illustrate our results, in fig. |^ we plot ^3(2;, y) at time t = for the choice a = 7 = 1 and m = 0.95. 



3 The Nonlinear Schrodinger Equation 

The NLSE with both attractive and repulsive nonlinearity has found many physical applications in 
several diverse areas including fiber optics, Bose-Einstein condensates and waveguides |Q. 

3.1 Case I: Attractive Nonlinearity 

The NLSE with attractive nonlinearity is given by {h = 2m = 1) 

iut + Uxx + u\u , (19) 
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Figure 2: V3{x,y,t = 0) vs. for m = 0.95, a = 7 = 1. 

where without any loss of generahty we have fixed the coefficient of the nonlinear term to be unity. 
As usual, one starts with the ansatz 

u{x, t) = r(^)e*(^(«)+"*) , ^ = x-vt, (20) 

which on substituting in eq. (|^) yields 



^ = ^(^^ + 4)' (21) 
2 



where prime denotes a derivative with respect to the argument and A, B are constants of integration. 
Thus the whole problem reduces to finding the solutions of eq. (p2|), after which 9 is easily obtained 



by using eq. (21) and performing one integration. 



The well known soliton solution of eq. (^) is 

r(^) = \/2 sech ^ , ^ = y - A = S = , = 4(n - 1) , (23) 

which is valid only for n > 1. It may be noted that a somewhat more general solution with arbitrary 
amplitude a is easily obtained, since if u(x,t) is a solution of the NLSE, then au{ax,a'^t) is also a 
solution of the same equation. 
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The two simplest, periodic, cnoidal traveling wave solutions of eq. (p^) are {^i = x — vit) 

ri(0 = V2dn ei , ^1 = ^ , A = , 5 = 4(1 - m) , vj = 4{n + m - 2) , (24) 

ri(^) = V2^cn 6 , Oi = ^ , ^ = , B = -4m(l - m) , = 4(n + 1 - 2m) . (25) 



In the limiting case m = 1, one recovers the familiar soliton solution (^31). 

We shall now show that suitable linear combinations of the solutions (123) and (1251) are also solutions 



of eq. (^). Consider first the solution (24). Our solutions consist of adding terms of the kind given in 
this equation but centered at p equally spaced points along the period 2K{m), where p is any integer. 
The p-point solution is 

rp{x, t) = V2Y,di, di = dn[^p + ^iinH^M, ^] , = {x - Vpt) . (26) 
1=1 P 

Clearly, p = 1 is the original solution, but for any other p, we have solutions of period 2K(m)/p. 

In order to verify that expression ( p^ ) is indeed a solution of eq. (p2|) , one needs the identities 

p p 



{Y^d.f = j2d^ + A{p,: 

i=l i=l 

ij2 d^f = jZdi + C{p, m) ^ + ^(^^ 



i=l i=l 
P P P 

■ . ^, 

i=l 1=1 i=l 

P P 

w? ^ SiCiSjCj = E{p, m) ^ (if + F{p, m) , (27) 

i<j 1=1 

which can be easily established by following the procedure discussed in ref. The general expression 
for the velocity Vp is 

v'^ = 4[n + m-2- C{p, m) - 2E{p, m)] . (28) 
Some explicitly computed values of the constants C{p,m) and E{p,m) are 

d>mq 



C{2, m) = 4E{2, m) = 4Vl - m , C(3, m) = AE{3, m 



C(4, m) = AE{A, m) = 4t(2 + t + 2p) , (29) 
where q is given by eq. (^) and i is given by 

t=(l-m)i/^. (30) 

On the other hand, for any p at m = 0, C{p, 0) = 4:E{p, 0) = '^^^ and at m = 1, C{p, 1) = E{p, 1) = 

^2 



0. It then follows from eq. ( P8|) that the solution as given by eq. ( ^6|) is valid only if n > 2p , and 
in this case v"^ changes from 4(n — 2p^) to 4(n — 1) as m goes from to 1. 



For odd p, using the solution (P5[), we obtain the following solution of the NLSE (|2^ ) by linear 
superposition 

/ , / r 4:(i-l)K(m) 1 , , 

rp[x,t) = \/2m2_^Ci , Cj = cn[r]p H ,m\ , rjp = [x - Vpt) . (31) 

i=i P 
In order to verify that (^) is indeed a solution to the NLSE eq. (22) one needs the identities 

p p 



i=l i=l 

{jZc^f=jZci + H{p,m)j2ci+I{p, 



m] 



i=l i=l i=l 

P P 

m- ^ SidiSjdj = J {p, m)^cf + K{p, 



.2\ ^-^ _ ^c^f + K{p,m) , (32) 

i<j 1=1 

which can be established following the procedure discussed in ref. Q. The general expression for the 
velocity Vp is 

Vp = 4[n + 1 — 2m — mH{p, m) — 2J{p, m)] . (33) 

Some explicitly computed values of the constants H{p, m) and J{p, m) are 

q + 2 



mH{3,m) = 4J(3,m) = -Aq[- 



(34) 



(l+g)2 l-g2J 

where q is given by eq. Thus v"^ varies from 4(n + 9) to 4(n — 1) as the elliptic modulus parameter 
m changes from to 1. 

On the other hand, for even p, we have obtained the following solution of the NLSE eq. (p^): 

p 



rp{x,t) = V2 ^[di- di+i] 

i odd 



(35) 



In order to verify that eq. (35) is a solution of the NLSE eq. ( |2^ ) one needs the identities 

P \ 2 P 



iZidi-d^+i]) =Y^dj + P{p, 

i odd 1=1 



4 P 



Y,[d,-d,+M =J2df + Lip,m)Y,d^ + M{p, 



m) 



m 



E 



SiCi SjCj 



E 



SiCi SjCj 



N{p,m)j2dl+Q{p, 



(36) 



i=l 



■ i+j even i+j odd 

which can be established following the procedure discussed in ref. Q. The general expression for the 
velocity Vp is 



n + m — 2 — L{p, m) — 2N{p, m) 



(37) 



Some explicitly computed values of the constants L{p, m) and N{p, m) are 

L(2, m) = AN{2, m) = -4^1 - m ; L{4,m) = AN{4,m) = -4i{2 - i + 2P) , (38) 



where i is as given by eq. (|30|). Thus for p = 2 [4], varies from 4(n + 4) [4(n + 16)] to 4(n — 1) 
as m changes from to 1. Generalizing the results in eqs. (^) or (|35|) , one finds that varies from 
4(n +p^) to 4(n — 1) as m changes from to 1. 

3.2 Case II: Repulsive Nonlinearity 

The NLSE with repulsive nonlinearity is given by 

iut + - u \ u\'^= , (39) 



We again start with the ansatz given by eq. (^0|) and on following the same steps as given in eqs. (20) 



to ( p2D it is easily seen that the 6 equation (eq. (21)) is the same as before while the r equation is 



almost the same except for the sign of the r'^ term. In particular, the r equation is now given by 



The well known soliton solution to this equation is 

r(^) = ^/2 tanh ^ , 6* = y , A = , B = -A , v"^ = 4{n + 2) . (41) 
The simplest, periodic, cnoidal traveling wave solution to eq. (^) is 

ri(C) = \/2msn , 6*1 = , A = , B = -4m , vj = 4:{n + I + m) . (42) 



For odd p, using the solution eq. (42), we obtain the following solutions of the NLSE eq. ( ^4 ) by 
linear superposition 

, , / r 4(i-l)K(m) 1 , X / X 

rp{x,t) = \/2m2_^Si , Si = sn[r]p H ,mj , rjp = [x - Vpt) . (43) 

i=i P 
In order to verify that eq. (43) is a solution to the NLSE eq. ( |40| ) one needs the identities 

(E5.)^ = E^1 + ^(p,^ 

i=l i=l 

(E ~Si)^ = jl~st + S{p, m) E s1 + T{p, ^ 



j=l i=l 
V V 

i=l i=l 1=1 
P P 



E CidiCjdj = U {p, m) E + Y{p, m) , (44) 

i<j i=l 



which can be estabhshed by fohowing S. The general expression for the velocity Vp is 



Vp = 4[n + 1 + m + mS{p, m) — 2U {p, m)] . 

Some explicitly computed values of the constants H{p, m) and J(p, m) are 

1 

mS{3,m) = — 4C/(3, m) = 4m[ 



m 



(45) 



(46) 



where q is given by eq. (^). Thus v'^ changes from 4(n + 9) to 4(n + 2) as m changes from to 1. 

For even integer p, the linear superposition of elementary solutions does not work. However, 
remarkably enough we find that the products of elementary solutions are also solutions. For example, 
the solution for p = 2 is 

r2{x,t) = V2msiS2 , (47) 

and the corresponding velocity is given by = 4(n + 4 — 2m). Generalization to higher even values 
of p is straightforward. 



3.3 Solutions With A^O 

It may be noted that since A = for all the solutions discussed so far, the expressions for 9 were 
rather trivial. One way of obtaining a solution with j4 7^ is to start with the ansatz 

r^iC) = 2dn^ ^ + a , (48) 



where a is a constant. It is easily checked that (48) is a solution to the NLSE (22) provided 



2 

a = -{n-2 + m- —) = Aa[a'^ - {n - —)a + 2(1 - m)] , 



B 



-3a^ + 4(n )a - 4(1 - m) 



(49) 



Starting from the solution (|48| ) we can obtain a class of solutions by an appropriate linear super 
position. For example, the 2-point solution is 



(50) 



It is easily checked that this is indeed a solution provided 

a 



2 

-{n-2 + m-—), B = -3a^ + 4(n - j)a - 16(1 - m) 



A' 



— in + 8(1 — m)a + 8(1 — m){n + 2 — m 



(51) 



Generalization to arbitrary p is straightforward. 
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4 The A(/)4 Model 

The kink (domain wall) solutions to the X(j)'^ field theory in (1 + 1) dimensions 

(pxx - (l>tt = A0((/>^ - a^) , 
have been widely discussed in the literature P]. The famous static kink solution is 



(x) = atanh(A/ A/2 ax 



from which the time-dependent solution 

(/)(rE, t) = atanh 



A 



■ a{x — vt) 



(52) 



(53) 



(54) 



2(1 -1-2) 

is immediately obtained by Lorentz boosting. Therefore, to begin with, we shall discuss only the static 
periodic kink solutions in case < 1 (and A > 1). Later, we shall discuss time dependent solutions 
with > 1 (or A < 1). 

4.1 Static Periodic Kink Solutions 



It is well known that the static periodic kink solution to the field eq. (52) is 

(l)i{x) -- 



2m 



■a sn(r/i, m) , rji 



A 



■ ax . 



(55) 



1 + m y 1 + m 

For any odd integer p, we find the following static kink solutions of the X(j)^ field theory by a 
specific linear superposition of the basic solution eq. (|55|): 



p 
1 

i=l 



2m aa ^ Si , p odd , 



(56) 



where Si,Ci,di are as given in eq. (14) with rjp = \/Aaax. In order to verify that eq. (^) is a static 
periodic kink solution to the X(j)^ theory field eq. (|5^), one needs the identity |^] 

p p p 



The constant a is given by 



2=1 



a 



i=l 



•\/l + m + 2mV{p, m) 
As an illustration, consider p = 3. V{2>, m) is given by 



mV{3, m) = 3 



m 



1-^2 
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(57) 



(58) 



(59) 




Figure 3: (j)i{x) (largest amplitude curve), (l)2{x) and </>3(x) (smallest amplitude curve) for m = 0.98, 
A = a = 1. 

where q is defined in eq. (M). Note that a varies from 1/3 to 1/^2 as m varies from to 1. 

For even integer p, the linear superposition of elementary solutions does not work. However, 
remarkably enough we find that even in this case, the products of elementary solutions are solutions. 
For example, the solution for p = 2 is 

4>2{x) = V2maasiS2 , (60) 

where si^2 are as defined in eq. (I) with e = aax. Using the identities derived in ref. it is 
easily shown that (pO]) is a static periodic kink solution to the field eq. ([5^) provided 



a = , / . (61) 
^2(2 - m) 

Note that a varies from 1/2 to 1/V2 as ni varies from to 1. Generalization to arbitrary even p is 
straightforward. It is thus clear that for arbitrary integer p, for static periodic kink solutions of Xcj)'^ 
theory, a will vary from l/p to l/\/2 as m varies from to 1. For the values m = 0.98 and A = o = 1, 
we plot (pi{x) from eq. (^), (p2{x) from eq. (|60|) and 03 (x) from eq. (p^) in fig. 0. 



4.2 Periodic Time-Dependent Kink Solutions 

While in the relativistic field theory context with A > 0, eqs. (^) and (|55| ) are the only solutions 
of A^^ field theory, in the condensed matter physics context, where velocity v can exceed velocity of 
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sound (optical modes), or for relativistic case with A < 0, one also has another soliton solution given 
by 

(j){x,t) = V2asech{l3{x-vt)) , [3 = a J , (62) 

y (f^ — Ij 

which is real if either A < (and < 1) or t;^ > 1 (and A > 0). 

The corresponding periodic soliton solutions to eq. ( |5^ are well known and are given by []To| 



<Ai (6 = y ^ « dn(6 , m) , 6 = a J _ J^^2 _ i) " > (63) 



l(^) = \/ o 1 Qcn(7?i,m) , rji = aJ — -^-^ -{x - vt) . (64) 

y 2m — 1 y (2m — l){vf — 1) 



Notice that the solution ( |6^ ) is valid only for 1/2 < m < 1 and both solutions (^) and ( plD are only 
valid for f ^ > 1 (and A > 0) or A < (and v"^ < 1). 

Appropriate linear superposition of solutions ( |63|) and (|6^ ) are also periodic time-dependent kink 
solutions. For example, by using the linear superposition of solutions (|63[), we have the following 
solution to the field eq. (^2|) , which is valid for any integer p: 

p 

cPpi^) = V2aaY,di, (65) 

i=l 



where di is as defined in eq. (^) with = a^J^^ja{x — Vpt). In order to prove that this is a solution, 
one needs the identity 

(E d^f = i24 + W[p, m)Y,d,, (66) 

i=l i i 

which is easily proved following ref. Q . Using this identity one finds that (^5|) is a solution provided 

" ~ [2-m + W{p,m)] ' ^^^"^ 



Some explicitly computed values of W{p, m) are 



W(2, m) = SVT^^ , 1^(3, m) = , WU, m) = 3i[2 + i+ 2p] , 

1 — q'' 

where i is as defined in eq. (|30|) and q is given by eq. (^). Further, W{p, 0) = — 1 , while W{p, 1) = 0. 
For any odd integer p, we also have the following solution to the field eq. (^2|) by linear superposition 

p 

4>{r]) = V 2m aa Cj , (69) 

2=1 
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where Cj is as defined in eq. (|T^ with r]p = ^•^_^ a{x — Vpt). In order to prove that this is a solution, 
one needs the identity 

d2c,f = j2cf + X{p,m)j2c^. (70) 

1=1 i 1=1 

Using this identity one finds that (|69|) is a solution provided 

2 _ \ /yix 

" ~ [2m-l + 2mX{p,m)] ' ^ ' 



For example, one can check 0] that 



6^2 



X(3,m) = -6(l-m + g) + — ^ , (72) 

with q being given by eq. (P), so that unlike the p = 1 case, this is an acceptable solution for all values 
of m (0 < m < 1). 

Similarly, for even integer p, we have solutions of the form 

p 

(j){C) = V2aaY,[di-di+i], (73) 

i odd 



where di is as defined in eq. (^) with = a^j^^^a{x — Vpt). Unfortunately, it appears that all these 
solutions are only valid in a very narrow range of values of m corresponding to real values of a. 

5 Sine-Gordon Field Theory 

In recent years, both sine-Gordon and Xcf)^ field theory have received considerable attention ^. 
In particular, sine-Gordon theory is the only relativistically invariant field theory having true soliton 
solutions. The equation under consideration is 

(t)xx - (ptt = sin (/) . (74) 

5.1 Static Soliton Solution by Linear Superposition 

The well known static one-soliton solution of this equation is given by 

0(x) =4tan-^e^'^ . (75) 

The corresponding time-dependent solution is easily obtained by Lorentz boosting and hence without 
any loss of generality we shall restrict our discussion to the static solution only (except when v'^ > 1). 
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The solution (75) can also be written in the alternative form 



sin(^Y^) = sech x . (76) 
The two corresponding periodic static soliton solutions are well known and are given by Q 

sin(^y^) =dn(x,m) , (77) 

sin(^y^) = cii[x / \/rn,m) , m > . (78) 

For any odd integer p, we obtain the following periodic static soliton solutions by linear superpo- 
sition: 

sin(^)=aX:^- (79) 
^ i=i 

where di is as defined in eq. (0) with r]p = ax, while a is given by 

«^ = T IT r TTr ^ • (80) 

[p + A[p, m) + mK[p, m)\ 



Here A{p, m) and i?(p, m) are as defined by eq. ( pTj ) and ( |44| ) respectively, and with this choice of a, 
one gets 

cos(^) = V^a^s^ . (81) 
^ i=i 

Note that use has been made of the identities |^ 

si[c2 + ... + Cp] + c.p. = , si[d2 + ••• + dp] + c.p. = , (ii[c2 + ... + Cp] + c.]5. = , (82) 



in proving that ( [79| ) is indeed a solution to the field eq. (74). For p = 3 the values of the constants are 



^(3, m) = 2q{q + 2) , mR{3, m) = 2{q^ - 1) , (83) 

where q is given by eq. (^), so that a = 1/(1 + 2q) changes from 1/3 to 1 as m varies from to 1. 
Another solution valid for any odd integer p is 

sin(M)=«^g^, (84) 
^ i=i 

where Cj is as defined in eq. (jlj) with r/p = axj ^pm. This solution is strictly valid only if m > 0. It is 
easily checked that this is indeed a solution to the field eq. (|7^) provided 

o? - ^ (85) 

\p + G{p, m) + R{p, m)] 
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Figure 4: sm[(j){x) /2] for p = 1, 2, 3 . Increasing p decreases the amplitude and increases the frequency. 



where G{p, m) and R{p, m) are as defined by eq. (^3|) and ( pOD respectively. Note that with this choice 
of a 



For p = 3, G(3, m) is given by 



cos(^) = af]5, . (86) 



where R{3,m) is as given by eq. (^), so that a = (m, > 0). 
Similarly, for p = 2 we have the solution 

sm(— — ) =a2_^di , 
^ i=\ 

where dj is as defined in eq. (|5|) with ^2 = olx. It is easily checked that this is indeed a solution of the 
field eq. (l?!) provided 

« = 1 ^ > (89) 

so that a varies between 1/2 and 1 as m changes from to 1. Note that with this choice of a 

0fx) 

cos( — ^ — ) = m,as\S2 ■ (90) 

In fig. ^, we have plotted sin[(^/2] vs. x for p = 1,2,3 corresponding to the right hand side of eqs. 
(|7^, (|8|) and (^). 

Another solution for p = 2 is 

sin(^) = a[di-d2] , (91) 
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where di is as defined in eq. d) with ^2 

= ax. It is easily checked that this is indeed a solution to the 
field eq. ( [7^ ) provided < m < 1 since a given by 

« = ^= , (92) 

1 — yl — m 

diverges at m = 0. 

5.2 Periodic Time-Dependent Solutions 

As in the Xcj)^ field theory case, in this case also we get solutions by linear superposition, which are 
only valid for > 1. In particular. Hp is any odd integer, the periodic time-dependent solution is 
given by 

sin(^) = V^af]5. , (93) 
^ i=i 



where Si is as defined in eq. (14) with rjp = a{x — vt)l\/v^ — 1. It is easily checked that this is indeed 
a solution of the field eq. ([7^) provided o? is again given by eq. (|80|) . Note that with this choice of a 

cos(^) = aX:'i.. (94) 
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Another solution for any odd integer p is 



sin(M)=«^5^, (95) 
^ i=i 

where is as defined in eq. (0) with r/p = a(x — vt) /\/m(f^ — 1). Thus, this solution is strictly valid 
only if m > 0. It is easily checked that this is indeed a solution of the field eq. (^) provided is as 
given by eq. (|85|). Note that with this choice of a, 

cos(^) = af:c,. (96) 
^ i=i 

Finally, for p = 2 we have the solution 

• /'^(^)^ /n'7^ 
sin( — —) = masiS2 , (97) 

where si^2 is as defined in eq. <^ with ,^2 = Ci{x — vt)/Vv'^ — 1. It is easily checked that this is indeed 
a solution of the field eq. ( |7^ provided satisfies eq. (p8|). Note that with this choice of a 

cos{^) = a j^d,. (98) 
^ i=i 
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6 Boussinesq Equation 

The Boussinesq equation is given by 

utt - + 3(n2)xx - Uxxxx = . (99) 
The periodic one soliton solution of this equation is known to be 

u{x, t) = -2a^dn2(a(x - vt)) + (ia^ , (100) 

where 

2(4 - 2m - 3/3)0^ = - 1 . (101) 

Thus v"^ <l \i (3 > 4/3 and > 1 if /? < 2/3. For 2/3 < /3 < 4/3, changes sign at some value of 
m (0 < m < 1). Note that in the limit m — > 1 and /? = 2, this solution goes over to the one-soliton 
solution 

u{x, t) = 2o? tanh^(Q(x - vt)) , = \J — ^ — , (102) 
Now consider the linear superposition 

p 

u{x,t) = -2a'^J2d'iia{x-vt))+(3a'^ , {p = 1,2,3...) . (103) 

i=l 

It is easy to check that this is an exact solution to the Boussinesq eq. ( p9| ) provided 

2[4 - 2m - 3/? - 6Ai {p, m)]a'^ = y"^ - 1 . (104) 

Here use has been made of the identity (^) with Ai{p, m) given by eqs. (^) and (p!o|). It may be noted 
that Ai{p,m) < 0. 

In this article we have shown that in view of the remarkable identities satisfied by Jacobi elliptic 
functions, a kind of linear superposition principle works for several nonlinear equations, some of which 
are associated with integrable systems, the others not. It would indeed be worthwhile to obtain such 
solutions for other nonlinear systems where elliptic functions play a role in the space of exact solutions. 
A question which comes to mind is how the solutions obtained in this paper are related to previously 
known solutions. At first sight, it would appear that our procedure has given new solutions, but a 
closer investigation reveals that our solutions are expressible in terms of previous solutions via a non- 
trivial generalization of Landen's formulas which connect Jacobi elliptic functions with two different 
modulus parameters [l^]. The reader is referred to Ref. ||ll[ for more details. 

One of us (A.K.) thanks the Department of Physics at the University of Illinois at Chicago for 
hospitality. We acknowledge grant support from the U.S. Department of Energy. 
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